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Figure  Captions 

Figure  1.  Axial  Pressure  Field  of  a 454  kHz  Signal  Radiated  by  a 3"  Diameter 
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Proiector  in  Fresh  Water  [r  = 1.5  yds,  a r = 7.35  x 10  Np, 

O 1 o 

B/A  = 4.9]. 

Figure  2a.  Axial  Pressure  Field  of  a 454  kHz  Signal  in  Fresh  Water. 

Figure  2b.  Axial  Pressure  Field  for  the  Second  Harmonic  of  a 454  kHz  Signal 
Generated  Via  Nonlinear  Self-Interaction  of  the  Fundamental  in 
Fresh  Water. 

Figure  3.  Finite-Amplitude  Absorption  Losses  Incurred  by  a 454  kHz  Signal 
in  Fresh  Water. 

Figure  4.  Far-Field  Finite-Amplitude  Absorption  Characteristics 

[20  Log10Oo  * SLo  + 20  Log10  </l  B k^c/) 

E SL*  + 20  LoglQ(2Tr/2  6 X 10-3/PqCo3)  = SLq  - ISO  dB 

in  water  (1  pbar  = reference  pressure) 

where  SL  = 20  Log(p  r / fl) 
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and  SL*  = 20  Log(p  r f / /l)  ; f in  kHz]. 

O 0 0 0 o 

Figure  5a.  Signal-Excess  Characteristics 

Figure  5b.  Up-Converted  Directivity  Characteristics 

Figure  6.  Approximate  Signal-Excess  Characteristics 
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Abstract 


•-‘■following  a review  of  the  literature  on  Parametric  Receiving  Arrays, 
the  problem  of  defining  their  performance  characteristics  under  saturated 
and  unsaturated  conditions  is  considered.  Basically,  the  problem  is  resolved 
by  establishing  equations  for  the  axial  field  of  a spatially  symmetric  pump 
wave  in  the  spectral  domain  via  Kuznetsov's  nonlinear  paraxial  wave  equation. 

As  a biproduct  of  this  analysis,  a simplification  of  terms  involving  the 
phase  of  pump  wave  in  these  equations  results,  upon  transformation  to  the 
time  domain,  in  a new  form  of  Burgers'  equation  for  a plane  piston  projector. 
Unlike  previous  forms  of  Burgers'  equation,  the  latter  combines  the  effect  of 
wave  interactions  in  the  near  and  far  field  regions  of  the  source.  Numerical 
comparisons  of  the  more  complete  spectral  equations  and  the  spectral  form  of 
the  new  Burgers'  equation  are  shown  to  be  in  good  agreement  with  experimental 
results  previously  reported  in  the  literature.  Approximate  solutions, of  these 
equations  are  also  derived.  The  three  methods  thus  established  for  represent- 
ing the  pump  field  are  then  used  to  derive  scaling  laws  for  parametric  receiving 
arrays,  which  clearly  show  the  limiting  effect  of  pump  wave  saturation  upon 
the  conversion  efficiency  of  the  up-conversion  process  as  the  pump  amplitude 


is  increased  indefinitely. 


-1- 


Introduc tion 

1.2 

In  his  work  on  the  scatter ing-of-sound-by-sound , Westervelt  ’ con- 
cluded that  two  overlapping  perfectly  collinated  plane  waves  of  finite- 
amplitude  would  only  give  rise  to  scattered  waves  when  propagating  -n  the 
same  direction  - a result  in  keeping  with  the  anharmonic  resonance  conditions 
discussed  by  Landau  and  Lifshitz.  This  work  led  naturally  to  his  funda- 
mental  paper  on  the  Parametric  Array  where  Westervelt  deduced  that  two 
perfectly  collinated  coterminous  progressive  f iniLe-amplitude  plane  waves 
of  angular  frequencies  and  traveling  in  the  same  direction  in  an 

unbounded  nondispersive  fluid  medium  would  interact  to  produce  highly 
directive  intermodulation  spectral  components;  that  of  frequency  , 

referred  to  as  the  "difference-frequency  signal"  being  the  lowest  in  the 
spectrum.  Since  he  was  primarily  interested  in  the  generation  of  low 
frequency  waves,  Westervelt ' s^  analysis  required  that  the  primary  frequencies 
be  very  nearly  equal.  This  requirement,  which  makes  possible  the  process 
of  frequency  "down-conversion,"  is  basic  to  parametric  transmitting  arrays, 
although  of  course  upper  sideband  components  are  also  formed. 

Alternatively,  when  viewed  analytically,  simultaneous  radiation  by  the 
source  of  frequencies  and  with  equal  finite-amplitudes  is 

equivalent  to  sinusoidal  modulation  of  a finite-amplitude  carrier  wave  or 
, Y,pump"  wave  of  frequency  (jJq  = ^y  a signal  of  frequency  . /2  , 

where  fi  = <d  - oj2  . It  follows  therefore,  from  the  inherent  "quadratic" 
nonlinearity  of  the  medium  that  the  spectral  components  of  the  modulating 
envelope  squared  and  in  particular  the  difference-frequency  ft  , will  be 
amplified  or  "pumped"  at  the  expense  of  the  carrier  until  the  amplitude  of 
the  latter  is  reduced  by  this  and  other  absorption  losses  to  a level  where 


-2- 


g ■ 


I 

. . 


•«» 

ton 


IL 


it  can  no  longer  sustain  a nonlinear  interaction.  The  reader  should  note 

at  this  point  that  the  interaction  is  not  confined  to  suppressed  carrier 
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modulation,  as  exemplified  by  the  work  of  Eller  and  Merklinger  on  more 
general  forms  of  modulation.  Viewed  in  this  manner  however,  it  is  clear 
that  the  process  of  parametric  amplification  can  be  interpreted  in  terms 
of  the  concept  of  finite-amplitude  self-demodulation  introduced  by  Berktay^ 

I 

and  justified  experimentally  by  Moffett,  Westervelt,  and  Beyer. 

Likewise,  in  the  converse  process  of  frequency  "up-conversion,"  para- 
metric receiving  arrays  are  formed  by  projecting  a finite-amplitude  pump 
wave  of  angular  frequency  u>  into  the  medium  to  serve  as  the  carrier  wave 
for  a weak  incoming  signal  of  frequency  ft  , where  to^/ft  >>  1 . In  this 
instance  the  pump  field  is  augmented  by  the  spatial  component 
of  the  signal  along  the  pump  axis,  assuming  of  course,  that  both  the  pump 
wave  and  this  component  are  traveling  in  the  same  direction.  The  resulting 
nonlinear  interaction  gives  rise  to  an  intermodulation  spectrum  as  before, 
the  sum  and  difference  frequency  components  + ft  being  of  greatest 

interest.  Since  u>  /ft  is  assumed  to  be  considerably  greater  than  unity 

o 

however,  these  sidebands  are  now  in  close  spectral  proximity  to  the  pump 
frequency,  but  unlike  the  latter  their  directivity  is 

equivalent  to  that  of  a virtual-end-fire  line  array  of  length  L/A^  (in 

wavelengths  of  the  signal  frequency)  where  L is  the  distance  from  the  pump 

projector  along  its  axis  at  which  a receiving  hydrophone  resonant  at  ojq  + ft 

or  oj  - ft  is  located.  Upon  reception  the  "up-converted"  signal  is  fed 
o 

to  a low  pass  filter  to  remove  the  pump  frequency  and  recover  the  signal  of 
frequency  ft  . 
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Although  implicit  in  Westervelt's  work,  ’ the  process  of  Parametric 

reception  was  identified  and  made  explicit  by  the  extensive  theoretical 

9 

and  experimental  investigations  of  Berktay  who  in  cooperation  with  Al- 

Ternim!10-12  considered  the  practical  implications  of  the  up-conversion 

13  1 A 

process.  Subsequent  experimental  work  ’ has  been  directed  to  long  wave- 

13 

length  up-conversion  in  fresh  water  lakes  and  to  the  consideration  of  arrays 

of  parametric  receivers,  thus  involving  significant  practical  extensions  of 

the  original  scaled  laboratory  experiments.  Further  theoretical  extensions  ’ 

have  also  been  made  to  provide  a more  precise  description  of  the  pump  fields 

radiated  by  practical  sources  and  the  resulting  effect  of  such  refinements 

upon  the  analytical  form  of  solutions  for  the  up-converted  fields.  More 

17  18 

recently  Goldsberry  and  McDonough  have  derived  optimum  operating  con- 
ditions for  parametric  receiving  arrays  from  systems  analyses  based  on 

Berktay  and  Al-Temimi's  analytical  model10  for  a spherically  spreading 

19 

pump  wave.  With  the  exception  of  a preliminary  study  by  Bartram  there 

has  been  no  systematic  study  of  the  effect  of  finite-amplitude  absorption 

on  the  performance  of  parametric  receivers,  which  although  insignificant 

at  low  pump  amplitudes,  ultimately  determines  the  maximum  achievable  efficiency 

of  these  arrays  when  the  pump  wave  becomes  saturated.  In  order  to  provide 

, 20 

a more  complete  analysis  of  this  effect  we  now  consider  how  Kuznetsov  s 
nonlinear  paraxial  wave  equation  (which  on  account  of  its  parabolic  form 
is  ideally  suited  to  numerical  solution)  can  be  used  to  define  the  acoustic 
field  of  a pump  wave  radiated  by  a plane  piston  projector  ’.nder  both 
saturated  and  unsaturated  conditions.  Following  this  analysis  we  will 
then  determine  the  effect  of  saturated  pump  waves  on  the  conversion 
efficiency  of  parametric  receiving  arrays,  and  hence  deduce  scaling  laws 
to  define  their  performance  capabilities. 
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1 . Kuznetsov  s Equation 

In  this  section  we  wish  to  investigate  the  distortion  of  progressive 

20 

finite-amplitude  waves  in  an  unbounded  nondispersive  fluid  via  Kuznetsov's 

nonlinear  paraxial  wave  equation,  which  it  should  be  noted,  is  a more  general 

form  of  the  inviscid  paraxial  wave  equation  previously  derived  by  Zabolotskaya 

and  Khokhlov.'  ’ Introducing  a rectilinear  Cartesian  coordinate  system  and 

considering  the  propagation  of  progressive  plane  finite-amplitude  waves 

20 

along  the  positive  z axis,  Kuznetsov's  equation  can  be  expressed  as, 


,2  ~ 

- (1/2  kn)VZ  u 
3t3z  ft  xy 


3t 


(B£okft)u  3t  + “ft  3t2  } 


wnere 


and 


v/vo  , 
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e = v /c 
o o o 


k = ft/ c , T = ft(t 
y o 


xy 


3x 


3y 


(1) 


-)  (2) 


(3) 


In  this  notation  v denotes  the  particle  velocity  resulting  from  propagation 

of  an  acoustic  disturbance  in  the  fluid,  v being  its  peak  value  at  the 

o 

source;  p is  the  static  density,  c is  the  small-signal  speed-of-sound , 
o o 

e is  the  peak  Mach  number  at  the  source,  and  6 is  the  second-order  non- 
o 

linear  coefficient  of  the  fluid;  ft  is  an,  as  yet  unspecified,  angular 
frequency,  and  is  the  corresponding  small-signal  thermo-viscous 

attenuation  coefficient. 

Assuming  that  an  axially  symmetric  disturbance  is  established  in  the  field 

via  radiation  from  a baffled  piston  source  of  radius  a , it  is  convenient  to 

normalize  the  variables  x and  y with  respect  to  a and  combine  them  in 

23 

a single  variable  £,  defined  by  Rudenko,  Soluyan,  and  Khokhlov  as, 
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£ = (x/a)2  + (y/a)2 


(A) 


Expressing  the  x and  y derivatives  in  terms  of  £ we  obtain, 


JL  !£  JL  JL  = 31  JL 

3x  3x  3£  ’ 3y  3y  ^ 


hence,  from  Equations  (3)  and  (4) 


xy 
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3 ' 
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3€. 

1 9y 

/ 

2 

2 1 
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laJ 

3C  P 3C  J 

(5) 


Substituting  Equation  (5)  in  Equation  (1)  thus  gives, 


- V2/2)_1  ik 


« |)  - £ It  + "B  0 } • (6) 


Normalizing  z with  respect  to  a reference  ’Rayleigh  distance'  rQ  kQa  /2 

where  k = u /c  , and  the,  as  yet  unspecified  frequency  w is  assumed  to 
o o o 

be  radiated  directly  by  tne  source,  Equation  (6)  becomes, 


a2 

3 u 

3x3R 
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3C 


_3_ 

3t 
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where  R - z/r  , r = k a2/2  , k - w /c  , - Se Q’ Vo  • <8) 


Alternatively,  in  terms  of  the  'stretched  coordinate  system'  introduced 
by  BlackstockV  Equation  (7)  becomes, 
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Acoustic  Reynolds  number. 


(10) 


It  is  clear  from  inspection  of  Equations  (7)  and  (9)  that  for  one  dimensional 

, J 24,25 

waves  they  reduce  to  the  plane  wave  form  of  Burgers  equation. 

Returning  to  Equation  (7)  and  expressing  u in  terms  of  its  Fourier 

transform  u we  obtain, 

co 
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3u 
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u(R,£,t)  = 


1 r°° 

— / u ( R , „ ) e 

2tt  co  w 


i (w/R)t 
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(12b) 


and 


u *u 

w w 


tt-  / u / u / d (w/R) 
27T  >-co  w-w  w 


(12c) 


a = an(w/R)‘ 
w R 


(12d) 


If  the  signal  of  angular  frequency  w^  radiated  by  the  source  is  sinusoidally 

modulated  by  a pure  tone  of  angular  frequency  ^ then  for  integral  values 

of  N = w /R  > 1 , w assumes  the  discrete  values  w = nR  and  Equation  (11) 
o — n 

becomes  the  infinitely  coupled  set  of  partial  differential  equations  given  by, 
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(13) 


where  u*  , denoting  the  complex  conjugate  of  u , appears  in  Equation  (13) 
m m 

* 

because  u(R,£,t)  is  a real  valued  function  and  consequently  u_m  = um 
from  Equation  (12a).  Obviously,  if  the  source  waveform  is  unmodulated,  we 
simply  set  N = 1 . 

At  this  point  we  note  that  since  the  left-hand-side  of  Equation  (13)  is 
Similar  to  the  time  dependent  form  of  Schrodinger’s  equation,  we  can  express 
the  axially  symmetric  spectral  amplitudes  un(R,£)  as  a Laguerre  polynomial 
expansion  given  by, 

(n/N)C/l  ~ iR  i ^k(R)L,  (?)  (14) 


u (R,£)  - e - r , 

n kc02  n k 


where  the  Laguerre  polynomials  L,  (£)  obey  the  following  relationships, 

iC 
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CL’k(C)  - kLk(U  = "kLk_l(0 


(15c) 


L (£)  = 1 ; LAO  = -£  + 1 ; L9(0  = £ - 4£  + 2 , etc.(15d) 

O 1 ^ 


Substituting  Equation  (14)  in  Equation  (13)  and  employing  Equations  (15a) 
through  (15c)  to  simplify  the  resulting  expressions,  we  obtain  after  some 


manipulation, 
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In  order  to  solve  Equation  (16),  we  require  the  initial  values  ^(0)  . These 
can  be  obtained  from  Equation  (14)  if  the  boundary  value  u (0,0  is  pre- 


scribed. Thus,  by  means  of  the  orthogonal  relationship  defined  by  Equation 


(15a)  we  have, 
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If,  for  example,  u^(0,O  = 1,  0 £ S £ 1 
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and  N = 1 then, 
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In  general,  the  tripple  sums  on  the  right-hand-side  of  Equal  . (16) 

are  very  difficult  to  evaluate,  particularly  as  Equations  (17a)  and  (17b) 
cannot  be  reduced  to  simple  closed  form  expressions  if  k,  s,  £ are  non 
zero.  For  any  value  of  k however,  since  the  largest  terms  in  the  double 
summation  over  s and  £ are  those  for  which  s = £ , we  can  neglect  the 
terms  s i £ over  small  increments  of  the  distortion  distance  O^FfR)  , 
where  F (R)  is  a function  which  has  yet  to  be  defined.  Moreover,  since 
the  field  is  most  intense  along  the  beam  axis,  the  particular  values 
s = £ = 0 are  dominant,  thus  giving  the  resulting  approximation  to  Equation 
(16),  valid  for  small  values  of  O^F(R)  , 
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oo 
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(21) 


(22a) 


1 + 


2m/N 

1 + r" 


(22b) 


Inspection  of  Equation  (21)  shows  that  the  zero  order  Laguerre  mode  (i.e. 
k = 0)  is  now  independent  of  the  higher  order  modes  (i.e.  k = 1,2,  . . .), 
which  can  all  be  derived  from  it  a-posteriori . Physically,  this  implies 
that  the  axial  field  is  predominantly  defined  by  the  zero  order  Laguerre 
mode,  which  in  turn  is  primarily  responsible  for  distortion  of  the  off-axis 
field,  the  latter  having  little  effect  on  the  axial  field  over  small  dis- 
tortion distances. 
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Furthermore,  since  small  increments  of  the  'scaled*  range  O^FCR)  imply 
large  increments  in  actual  range  R for  <<  1 , the  model  ensures  that 
all  Laguerre  modes  other  than  zero  will  be  insignificantly  small  under 
unsaturated  conditions.  Thus,  for  a monofrequency  source  (i.e.  N - 1)  the 
unsaturated  harmonic  directivity  functions  formed  by  self  interaction  of 
the  fundamental  frequency  component  in  the  medium  are  given  by  Equation  (14) 
as  Gaussian  beam  patterns, 

e-n£(l  - iR)/(l  + R2)  j | £-nC(l  ~ iR)/(l  + R2)j 

= e“2n^/(l  + R2)  (23a) 

- e‘2n^/R2  = e~  f(koa  Sin  9)2  , R » 1 . (23b) 


D (R.O  = 
n 


Thus,  in  keeping  with  Rudenko,  Soluyan  and  Khokhlov's  investigation,  the 
model  is  capable  of  defining  the  major  lobe,  but  not  the  minor  lobe  structure 
of  the  actual  field.  Under  saturated  conditions  as  more  and  more  higher- 
order  Laguerre  modes  are  included  in  Equation  (14) , the  harmonic  directivity 
functions  become, 


Dn(R.O 


u (R,£)u  (R,£) 

n n 

u (R,0)u  (R,0) 
n n 


1 + 


OO  75 

Z Z <|;k(R)>r  (R)L.  (OL AO 

k.fcd  n n k & 

oo  k 

Z Z /(R)/  (R)L.  (0)L.(0) 
k,^l  n n k x, 


-2n£/(l  + R ) 


(24) 
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Under  unsaturated  conditions,  ^ -*■  0 for  k i 0 , so  that  Equation  (24) 

reduces  to  the  form  of  Equation  (23)  as  required.  With  increasing  values 

of  0^  therefore,  as  the  number  of  terms  in  the  summations  of  Equation  (24) 

increase,  the  directivity  functions  should  broaden  in  accordance  with  the 

27 

experimental  results  of  Shooter,  Muir  and  Blackstock. 

Since  the  zero-order  Laguerre  mode  can  be  evaluated  independently  of 
the  higher  order  modes  for  the  model  under  consideration,  it  is  expedient 
to  consider  the  differential  equations  for  this  mode  alone.  These  are 
given  by  Equation  (21)  as 
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(25) 


In  Appendix  A it  is  shown  that  under  unsaturated  conditions  the  second 

harmonic  field  in  an  inviscid  fluid  obtained  from  Equation  (25)  by  the 

method  of  successive  approximations  is  identical  as  a, r ->0  to  the 

1 o 

approximate  solution  previously  derived  from  the  inviscid  form  of 

23 

Equation  ( 9 ) by  Rudenko,  Soluyan  and  Khokhlov. 

Since  ip is  complex  it  is  convenient  to  express  it  in  terms  of  its 
real  and  imaginary  parts  ^ and  respectively  where, 


o 

n 


TO 

*n 


(26) 


Substituting  Equation  (26)  in  Equation  (25)  and  equating  the  real  and 
imaginary  terms  respectively  we  thus  obtain, 
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(27) 


With  the  exception  of  the  "diffraction-induced  spreading  losses"  on  the 

2m/N  1 


1/1  + 


in 


left-hand-side  of  Equation  (27)  and  the  coefficient 

1 + R‘ 

the  summation  on  the  right-hand-side,  the  coupled  harmonic  modes  are  similar 

28  . 

to  those  previously  considered  by  Bellman,  Azen  and  Richardson  in  their 

numerical  analysis  of  the  plane  wave  form  of  Burgers'  equation.  Solving 

these  equations  for  a monofrequency  source  (i.e.  N = 1)  by  a predictor- 

29 

corrector  method  of  the  Adams  Moulton  type,  we  examined,  among  other 

27 

cases,  the  experiment  conducted  by  Shooter,  Muir  and  Blackstock.  In  this 

instance  a 3"  diameter  plane  piston  projector  operating  at  A5A  kHz  in  fresh 

water  at  17.8°F  was  driven  at  source  levels  from  100-135  dB  re  1 Ubar  at  1 yd. 

27 

Assuming  the  same  parameters  as  those  previously  chosen,  we  have 

o^/f 2 = 2.6  x 10_lANp/m  at  17.8°F,  rQ  = 1.5  yds,  and  B/A  = A. 9 ; 

consequently,  we  obtain  aj^r0  = 7.35  x 10  Np.  Varying  in  accordance 

with  the  range  of  source  levels  mentioned  above,  wa  computed  ip  at  par- 

27 


ticular  distances  from  the  source  selected  for  the  experiment,  giving  the 
axial  pressure  functions  for  the  fundamental  frequency  component  depicted 
in  Figure  1.  It  can  be  seen  that  our  results  are  in  excellent  agreement  with 


measurements  except  at  ranges  close  to  the  source  (i.e.  0.76  yd  and  1.8  yd) 

27 

where  an  examination  of  Figure  6 in  Shooter,  Muir  end  Blackstock's  article 

shows  that  the  experimental  curves  at  these  distances  are  not  consistently 

scaled  (as  they  are  at  all  other  ranges)  to  give  an  extrapolated  peak  source 

level  at  1 metre.  For  comparison,  computed  values  of  the  fundamental,  and 

second  harmonic  are  shovn  in  Figures  2a  2b,  and  in  Figure  3 the  "extra 

decibel  loss"  is  given  as  a function  of  R at  different  source  levels 

corresponding  to  those  of  the  experiment.  Although  the  harmonic  directivity 

functions  for  this  example  are  not  shown  here,  they  can  be  computed  via 

Equations  (21)  and  (24)  from  the  data  already  obtained.  Our  results  have 

clearly  shown  however,  chat  the  axial  field  is  accurately  represented  by 

ij>°  (R)  , thus  justifying  the  assumptions  made  in  deriving  Equation  (21). 
n 

Since  the  real  and  imaginary  parts  of  Equation  (27)  are  coupled  on  the 

2 

left-hand-side  by  terms  whose  coefficient  (1/1  + R ) approaches  zero  as  1/R 
for  R >>  1 , and  since  the  coefficient  [ 1/1  + - 2m;/N  on  the  right-hand-side 


1 + R 


approaches  1 as  R increases  indefinitely,  we  now  propose  to  neglect  these 
terms  so  that  recombining  the  veal  and  imaginary  parts  of  we  obtain  an 

approximate  form  of  Equation  (25)  given  by, 


dUJ  wa  \ n-1  00 

—2.  + a r + ip°  = i - — V + 2 Z • (28) 

dR  ^ ° ! + r2  " IwJiAj  m m=1  n-hn  m 


In  this  equation  the  effect  of  diffraction  on  the  acoustic  field  is  maintained 

2 

by  the  spreading-loss  term,  whose  coefficient  (R/l  + R ) approaches  zero  as 

R approaches  zero  and  decreases  as  1/R  for  R >>  1 . Using  Equation  (28) 

to  rederive  the  unsaturated  solution  for  the  second  harmonic  formed  in  an 

i 23 

inviscid  fluid,  and  comparing  it  with  Rudenko,  Soluyan,  and  Khokhlov's 
solution  of  the  inviscid  form  of  Equation  (25),  as  given  by  Equation  (A  ) 
of  Appendix  A,  we  find  that  neglecting  phase  differences,  the  absolute  value 


>f  the  ratio  of  the  two  solutions,  designated  Q , can  be  expressed  as, 
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1 / , 2 _1 
2 / { Ln  ( 1 + R ) } + 4 { tan  _R }_ 


sinli  ^ R 


Evaluating  the  ratio  Q as  a function  of  R , it  can  be  seen  from  Table  I 


that  the  solution  of  Equation  (28)  for  the  second  harmonic,  approaches  but 


never  exceeds  10%  of  that  derived  from  Equation  (25) * It  can  also  be  seen 


from  Table  I or  from  inspection  of  Equation  (29)  that  the  difference  between 


the  two  solutions  tends  to  zero  as  R approaches  zero  or  infinity.  It  seems 


reasonable  to  assume  therefore,  that  Equation  (28)  will  provide  a good 


approximation  to  the  absolute  value  of  ^ . Further  confirmation  of  this 


conclusion  is  provided  by  Figures  1 and  2 for  the  case  of  Shooter,  Muir, 


and  Blackstock's  experiment,  where  values  of  the  fundamental,  and  second, 


harmonics  computed  from  Equation  (28)  are  seen  to  be  in  good  agreement 


with  the  results  obtained  from  Equation  (25)  based  on  the  same  numerical 


procedure. 


Using  an  approach  previously  adopted  by  the  author  in  applying 


Merklinger's  heuristic  plane  wave  analysis  to  the  spherical  wave  form  of 


Burgers'  equation,  it  can  be  shown  that  the  fundamental  frequency  component 
of  Equation  (28)  for  a monofrequency  source  (i.e.  N = 1)  may  be  represented 


by  means  of  an  approximate  expression  which  differs  slightly  from  a similar 


approximation  obtained  heuristically  by  Merklinger , Mellen,  and  Moffett. 


This  approximation  is  given  by, 
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33 

where  H and  Y are  zero  order  Struve  and  Neumann  functions,  respectively, 
o o 

As  shown  in  Figure  1,  results  obtained  from  Equation  (30a)  are  also  in  good 
agreement  with  those  computed  directly  from  Equation  (28),  thus  lending 
credibility  to  the  approximation. 

Another  interesting  result  inherent  in  the  form  of  Equation  (30)  is 
that  the  finite  amplitude  absorption  loss  incurred  by  the  fundamental  in 
propagating  through  the  medium,  designated  the  "extra  decibel  loss  by 
Blackstock3^  and  hence  denoted  EXDB  , can  be  expressed  as, 


Ola) 


(31b) 

Denoting  the  "extra  decibel  loss"  defined  by  Equation  (31b)  as  EXDB^  to 

indicate  that  for  fixed  values  of  a it  represents  the  maximum  finite-amplitude 

o 

absorption  loss  incurred  by  the  fundamental,  the  curves  of  Figure  3 were  computed 

showing  the  variation  of  this  parameter  with  Oq  and  a r . Since  0q  is 
related  to  the  source  level,  as  shown  in  the  caption  of  Figure  4,  these 
characteristics  can  be  used  to  obtain  the  "saturation  threshold"  of  a 
monofrequency  source  (i.e.  the  difference  between  a desired  source  level  and 

the  appropriate  value  of  EXDB^)  once  c^r  and  aQ  are  specified. 

Having  provided  some  degree  of  justification  for  Equation  (28),  we 

now  take  its  inverse  Fourier  transform  to  obtain  a new  form  of  Burgers' 


EXDB  = - 20  Log 
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equation  for  a plane  piston  source,  which  can  be  expressed  as 


9i 

3R 


1 + R 


- (a  /N)  ip  - (a  r ) 
O oT  o O 


i!|  , o 

9t 


(32) 


00 

where  ^(R,T)  = E 4>n(R)  e 

n=-ao 


(33) 


It  can  be  seen  by  inspection  that  if  R « 1 , Equation  (32)  approaches 

the  plane  wave  form  of  Burgers'  equation.  Alternatively,  if  R >>  1 , it 

assumes  the  spherical  wave  form  of  Burgers'  equation,  as  required. 

Assuming  further  that  the  axial  field  is  determined  primarily  by 

^°(R)  it  follows  from  Equation  (14)  that 
n 


u(R,o  ,t)  ~ iKR.T) 


(34) 


In  order  to  obtain  the  Fubini, 35,36  Fay,37-^0  and  asymptotic  far  field3^*41 
solutions  of  Equation  (32),  it  is  convenient  to  reexpress  it  in  terms  of  the 

9 / 

'stretched  coordinate  system  introduced  by  Blackstock.  This  transformation 
is  carried  out  in  Appendix  B yielding  the  scaled  form  of  Equation  (32)  pre- 
scribed by  Equation  (B7)  as, 

B-  <i/n>  hB-  {C  cosh<°/°0) } ■ ° <35> 

where  W = ip/ 1 + R^  (36) 

and  o/o  = F(R)  = sinh  ^R  . (37) 

o 

The  established  solutions35-^1  of  Equation  (32)  for  a monofrequency  source 
are  outlined  in  Appendix  B,  but  before  concluding  it  should  be  noted  that 
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the  unspecified  "distortion  function"  F (R)  previously  mentioned  in  the 
text  has  now  been  equated  to  sinh  "^R  in  Equation  (37).  Strictly  speaking, 
a more  precise  form  of  F(R)  for  the  case  of  Equation  (21)  would  be  the 
numerator  of  Equation  (29),  but  as  we  have  shown  in  Table  I,  this  differs 
only  marginally  from  sinh  ^R  . 

Having  tin  3 completed  our  discussion  of  the  field  equations  governing 
the  propagation  of  finite-amplitude  waves  radiated  by  a plane  piston  pro- 
jector, we  will  now  consider  how  these  equations  can  be  used  to  provide 
scaling  laws  for  parametric  receiving  arrays. 
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2 • Par a me  trie  Receiving  Arrays 

Since  the  analytical  model  that  has  been  developed  can  be  used  to 
define  the  field  of  a finite-amplitude  pump  wave  radiated  by  a plane  piston 
projector  under  saturated  or  unsaturated  conditions,  we  now  wish  to  consider 
the  up-converted  fields  generated  via  interaction  between  the  pump  and  a 
plane  wave  of  angular  frequency  ft  whose  wave  normal  intersects  the  pump 
axis  (i.e.  the  positive  z axis)  at  an  angle  0 . We  begin  by  expressing 
the  component  of  the  signal  along  the  pump  axis  in  ihe  form, 


u (z,t)  = Re 

s 


f i(ftt  - Xnz  cos  0)1 
u e ii 

s 

o 


- Re{us  ei(S!t  * V + 2V} 


- Re{us  e-V  + 1(flt  ‘ kfiz  + 2V>} 


(38) 


where  Xfi  = kQ  - iafi  , and  = X^  sin  (0/2) 


(39) 


We  now  let  t = ft ( t - z/c  ) = (ftt  - k z)  , R = z/r  , and  r = k a /2  so 

o it  o u U 

that  Equation  (38)  becomes, 
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If  oj  is  the  pump  frequency,  we  arbitrarily  chose  its  time  waveform  at 
o 

the  source  tc  be, 
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The  spectral  amplitude  of  the  component  of  the  signal  along  the  pump  axis  can 
thus  be  expressed  as 


*°(R) 


,c  -(a  r )R  + i(2Mnr  )R  . 
= U)  e 1 o 12  o ; 

s 

o 


hence  from  Equation  (26)  we  have, 

(43a) 
(43b) 

Likewise,  the  spectral  amplitudes  of  thp  pump  wave  at  the  source  can  be 
arbitrarily  expressed  as, 


li^(R)  = i{ )°  e~(0tiro)R  cos(2MnrQR) 
o 

<Jj°(R)  = "^s  e~(Ctir°)R  sin(2MftroR)  • 
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*N<°) 

= ^(0) 

(44) 

giving 

<(°) 

= 0 

(45a) 

and 

*N(°> 

= 1 . 

(45b) 

With  the  boundary  conditions  for  the  pump  field  prescribed  by  Equations  (45a) 
and  (45b),  we  can  now  solve  Equation  (27)  numerically  to  determine  the  up- 
converted  frequency  components  using  Equations  (43a)  and  (43b)  to 

A 

define  i|j°(R)  and  lp°(R)  throughout  the  field.  This  calculation  can  be 

repeated  using  different  values  of  N , (aNr0)  » and  a0  in  order  to  express 

the  'signal  excess'  ^n+1^s  aS  3 ^unction  t^ie  PumP  amplitude  via  the 

— o 

parameter  0 for  particular  values  of  (aj^r0^  and  R • These  calculations 


• • 
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must  then  be  repeated  for  different  values  of  0 in  order  to  determine  the 
spatial  directivities  of  the  up-converted  fields.  Sin^e  this  is  a very 
lengthy  process  it  is  expedient  to  take  advantage  of  the  fact  that  in  most 
applications  of  interest  the  signal  is  considered  to  be  very  weak  relative 
to  the  pump.  The  equations  for  the  up-converted  spectral  amplitudes  ^°+1 
can  thus  be  considerably  simplified  by  neglecting  terms  on  the  right-hand- 
side  of  Equation  (25)  other  than  those  that  involve  direct  products  of  the 
pump  and  signal  components  giving, 
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s. 


Under  unsaturated  conditions  the  pump  field  is  given  approximately  by  Equation  (25) 


as 
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e-(Vo)R 

^N(R)  “ I - iR 


(47) 
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Substituting  liquations  (42)  and  (47)  in  the  right-hand-side  of  Equation  (46a) 
we  thus  obtain, 
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Since  ip°  , (0)  = 0 , the  solution  of  Equation  (48)  enables  the  'signal  excess' 
N+l 

ip°  /{ h°  for  the  upper  sideband  to  be  expressed  as, 
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By  means  of  the  plane  wave  impedance  relationship,  we  can  reexpress  the  signal 

excess  as  ip°  /ip°  - p'/p'  so  that  in  terms  of  more  conventional  notation 

N+ Is  i s 

o +o 

Equation  (49)  becomes, 
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+ (Vo)R  sin  (0/2)}. (50) 


Inspection  of  Equation  (50)  shows  that  if  R « 1 it  approaches  Berktay  and 
Al.Temimi's  solution'lor  a plane  pump  wave,  and  likewise  if  R >>  1 it 
approaches  the  form  of  the  corresponding  spherical  pump  wave  solution, 
utilized  by  Goldsberry  ' and  McDonough  in  their  respective  systems  studies. 
Since  Equation  (50)  is  perfectly  general  and  can  easily  be  incorporated  in 
such  studies  it  should  prove  useful  in  cases  where  the  pump  projector  Rayleigh 


distance  is  large. 


7 
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Returning  now  to  Equation  (46b)  and  substituting  Equations  (42)  and  (47) 
in  the  right  hand  side  we  obtain. 


, ° 
dR 


[ aN-iro  " 1 - iR  J rN-l 


1 Ur 


= i 


N - 1 


( a 


l 2 J 's 


1/1  - iR 


1 + 


2/N 

1+R2 


1/1  + iR 
l+R 


- (a. , r + a r )R  - i(2MQr  )R 
; N o 1 o i!  o 


N - 1 


Ip~  U + 

s 

o 


1 - iR 
1 + iR 


i + R 
3 + R' 


- (a  r + a r )R  - i(2M0r  )R 
; N o 1 o o 

1 - iR 


~ l 


f 

2o 

N - 1 

o 

i « 

l 3 

-(a  r + anr  )R  - i (2Mfr  )R 
o e N o 1 o >2  o 


1 - iR 


(51) 


Hence  with 


p'/p'  the  solution  of  Equation  (51)  gives, 
s 

o 


p'/p' 


s 

o 


CO 

2o 

o 

R 

2 

sin  (Xnr  )R  sin  (0/2) 

H O 

CO 

l °J 

3 

/l  + R2 

(Xnr  )R  sin2 (0/2) 
U o 

-(a  r )R  + i{tan  ^R 
e - o 

-(Xnr  )R  sin2 (0/2) } . 
S2  o 


(52) 


Comparison  of  Equations  (50)  and  (52)  shows  that  with  the  exception  of  a slight 

• to  • 


change  of  phase  the  two  solutions  agree  to  within  a factor 


to 


0). 


co , 


-(a  r - a.r  )R 
e - o + o 


1 . 


Since  we  are  primarily  concerned  with  the  effect  of  a saturated 


pump  field  upon  the  up-converted  frequency  components  we  now  repeat  the  previous 
derivation  of  using  Equation  (30a)  in  place  of  Equation  (47)  to  detine 


the  pump  field.  In  this  manner,  neglecting  the  small  difference  between  , 
and  t|>°  , the  signal  excess  p'/p'  can  be  expressed  as. 


p 1 /p 1 ~ i 0 

' oj  2 

— o ° J 


a).  o I -(a.r  ) R + i tan  R rR 
+ o , e _ o 


/ 1 + r' 


<2Vo)R  dR1 


o rR’  -(2a  r )R" 

1+  -f  C -c 

^ ° / 1 + R”2 


’ oj . a -(a+r  ) R + i tan  1R  00  ? r1 

_ -j  _=  — — — — £ exp{i  m sin  (sin  0/2)} 

(jJ  2 O •'O 

o 


/ 1 + R 


J (2Xnr  R’)  dR’ 
m Si  o 


a rR’  -(2a  r )R" 
o e No 

2 J r 

° / 1 + R"2 


On  axis  this  becomes, 


p;/p;  - 

- o 0=0 


0J+  o 

Z.  o 


e-(a+rQ)R  + i tan  R rR 


/ 1 + r' 


fa  rR  e-(2Vo)R 

1+  -f  — =r  * 

/ ° / 1 + r"2 


Likewise,  if  D+(0)  denotes  the  spatial  directivity  functions  of  the  up- 
converted  fields  then  from  Equation  (53), 


E expfi  m sin  ^(sin2  0/2)) 


J (2Xor  R’)  dR’ 
m 21  o 

°o  (R'e-(2Vo)R" 


o . 2 

/ 1 + R" 


« « 
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Evaluating  Equations  (54)  and  (55)  numerically  we  obtained  the  results  depicted 

in  Figures  5a  and  5b  which  clearly  show  in  terms  of  'scaled'  coordinates  the 

effect  of  pump  wave  saturation  upon  the  up- conversion  process. 

Finally,  in  order  to  obtain  simple  numerical  estimates  for  values  of 

a r = a r <0.1,  we  make  use  of  weak  shock  theory  by  substituting 
N o o o 

i(;°(R)  z.  e”  aNro^/  [1  + Ln(l  - iR)  ] in  Equations  (46a)  and  (46b)  to  obtain, 


* * s 
+ o 


± 1 - i -f  R - Ei 

°o  I 2 > 


' 2 ~ 

o 

o 


- (a . r )R  + i tan  R - 2/o 
a Z o o 

/ 1 + R2 


0).  Ln  [1  - l -r-  R]  , \ . . -1_ 

Z 2 - a.r  ) R + i tan  R . . 

— e + o , o > 1 

to  — 7T  — o 

° / 1 + R2 


(57a) 


(i) . o 


->  i 


— e"(Vo)R  + 1 tan"  R , o R/2  < 1 


/ 1 + R^ 


(57b) 


[ O ] 2 1 2 r OqR  2 

Otp]  1 Lnl+fR  +4]tan"  -f-  -(a.r  )R  + i tan_1R  , 

— a — 1 —j e + o 

“oj  2/  1 + R2 

o R/2  > 1 (57 c) 

o 


Having  thus  obtained  approximate  expressions  for  the  signal  excess,  the 


"pump  excess"  p{/p_|.  can  likewise  be  approximated  by  means  of  Equations 
(3^a)  and  (50)  or  (57c). 

Since  the  "signal  excess"  as  defined  by  Equation  (57c)  is  virtually 

independent  of  to  /to  for  to  >>  ft  it  follows  that  "scaled  characteristics" 
+ o os 

giving  the  functional  dependence  of  (p_|/p'  ) exp{(a+rQ)R}  on  Oq  can 

— o — 

readily  be  constructed  for  different  values  of  the  "scaled  range"  R , as 
shown  in  Figure  6,  where  the  particular  choice  of  R equal  to  10  and  100 


respectively,  was  quite  arbitrary.  Using  such  characteristics,  the  signal 
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excess  can  readily  be  evaluated.  For  example,  in  the  case  of  a 100  kHz 

pump  wave  radiated  by  a 0.3m  diameter  projector  in  fresh  water,  with 

r = A.  7m  ct  = 2.6  x 10”4  Np/m  , and  a r = 1 . 23  x 10  3 Np  , Figure  6 
o 1 10 

gives  the  signal  excess  measured  by  a hydrophone  located  on  the  pump  axis 
at  R = 100  (i.e.  r = 470m)  for  0.01  £ aQ  £ 1 as, 


a 

Radiated  Pump  Power 

Signal  1 

o 

Watts 

dB 

0.01 

2 x 10"3 

-67 

0.1 

2 x 10-1 

-53 

1.0 

2 x 103 

-47 

It  is  clear  from  this  example  therefore,  that  increasing  the  radiated  pump 

power  to  enhance  the  signal  excess  reaches  a point  of  diminishing  returns 

when  a = 1 ; a conclusion  which  is  further  confirmed  by  examining  the 
o 

degradation  of  the  directivity  functions  as  increases. 
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Conclusions 


9 n 

Using  Kuznetsov's  nonlinear  paraxial  wave  equation  to  derive  coupled 

spectral  equations  for  the  axial  field  of  a inite-amplitude  pump  wave  radiated 

by  a plane  piston  projector,  we  have  shown  both  by  numerical  analysis  and  by 

successive  approximation  that  the  solutions  of  these  equations  are  in  good 

27 

agreement  with  well  established  experimental  results  under  saturated  and 
unsaturated  conditions.  We  have  also  shown  that  when  the  pump  field  is 
modulated  by  the  spatial  component  of  a weak  field  transmitted  along  the 
pump  axis  that  the  resulting  up-converted  frequency  components  can  be  amplified 
by  increasing  the  amplitude  of  the  pump  wave  until  a point  of  diminishing 
returns  is  reached  when  the  latter  exceeds  its  saturation  threshold.  In 
order  to  supplement  our  numerical  analysis  of  the  up-conversion  process,  we 
have  derived  analytical  expressions  for  the  signal  excess  and  pump  excess 
whicii  are  in  reasonably  good  agreement  with  our  numerical  results  for  values 
of  the  small-signal-pump-wave-absorption-loss  a0r0  less  than  0.1  np.  These 
closed  form  approximations,  which  reduce  at  low  pump  amplitudes  to  Berktay 
and  Al.  Temimi's10  spherical  pump  wave  solution,  are  sufficiently  simple  to 
be  readily  included  in  system  simulation  models  where  they  can  define  upper 
bounds  for  the  conversion  efficiency  and  directivity  of  Parametric  Receiving 
Arrays  as  the  pump  wave  undergoes  saturation. 
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Figure  2b.  Axial  Pressure  Field  for  the  Second  Harmonic  of  a 454  kHz  Signal 

Generated  Via  Nonlinear  Self-Interaction  of  the  Fundamental  in  Fresh 
Water. 


Figure  3.  Finite-Amplitude  Absorption  Losses  Incurred  by  a 454  kHz  Signal  in  Fresh  Water. 


Figure  5b.  Up-Converted  Directivity  Characteristics 


Figure  6.  Approximate  Signal-Excess  Characteristics 


APPENDIX  A 


For  an  unsaturated  monofrequency  source  of  finite-amplitude  (i.e. 

O < 1 , N = 1)  the  spectral  amplitude  of  the  fundamental  frequency 
o 

component  is  obtained  throughout  the  field,  to  a first  approximation,  by 
suppressing  the  right-hand-side  of  Equation  (22)  so  that, 


Vo  1 - iR  ^1 


-(a  r )R 
1 o 


<P°(0) 


giving  (R)  = "YTTr  * ( 

Substituting  Equation  (A2)  in  the  right-hand-side  of  Equation  (22),  the 
spectral  amplitude  of  the  second  harmonic  field  generated  by  self-interaction 
of  the  fundamental  frequency  component  in  the  medium  is  given  to  second-order 
by  the  equation, 


r2  i ,o 

dR  + a2ro  " 1 - iR  r2 


. -2a. r R 

ia  1 o 

2 — • (y  = 4a 

9 2 ’ 2 1 

2 [1  - iR] 


hence , 


^(R) 


-(a  r )R  D (a  r - 2a  r )R' 
ia  2 o cR  2 o 1 o 

o e *2 

2 1 - iR  1 - iR' 

'o 


o - l(2Vo>  ,1  - iR  1<2alro)y 


1 - iR 


-(4a  r ) R - i(2a  r ) 
0 1 o 1 o 

o e 

2 1 - iR 


| Eiti(2a1ro)(l  - iR)}  - E^i^a^)] 


where  Ei(x)  = 

X 

/— 00 

ey 

— dy  . 

y 

From  the  asymptotic 

limits  of  Ei(x) 

we  obtain,  for 

^(R)  * " 

r 

0 

' -(2a, r )R 
1 o 
e 

-(4a  r )R  ] 
1 0 
e 

4’ 

(1  - iR)2 

(1  - iR) 

(A5) 


f -(2a^r  ) R + 2i  tan 

t— 7771* — 


-(4a, r )R  + i tan  ^R ' 
1 o 


/ 


1 + R 


(A6) 


Alternatively,  as  a^rQ  -*  0 , 


_o  Ln(l  - iR) 
2 1 - iR 


!o  / {Ln(l  + R2)}2  + 4 {tan~1R}2 


i tan  1r 


- tan 


•l[  2 tan~1R  ] 
kn(l  + R2) 


1 + R 


(A7) 


This  inviscid  solution  is  identical  to  that  obtained  by  Rudenko,  Soluyan,  and 


Khokhlov . 


APPENDIX  B 


Combined  Plane-Spherieal-Wave  Burgers1  Equation  in  Stretched  Coordinates 

24 

Following  Blackstock  we  chose  to  redefine  y in  terms  of  a new 
variable  W , defined  as, 


W = \p/ 1 + R2  . 


(Bl) 


Substituting  Equation  (Bl)  in  Equation  (30) , we  obtain. 


- fa  /n/  1 + R2 
3R  o 


W |S  - (a  r ) & - 0 . 

9t  o o gT2 


(B2) 


f)  / 

Again,  following  Blackstock,  we  introduce  the  scaled  distortion  distance 

O = F(R)  in  Equation  (B2) , where  F(R)  has  yet  to  be  defined,  giving 
o 


3W 

do 


1/NF’ (R) 


/ 1 + R2  ] 


w U - [1/r  F'  (R)  ] 

3t  o 9t2 

r 


= 0 ; 


o /a  r 
o o o 


(B3) 


The  function  F(R)  can  now  be  defined  by  setting  F'(R) 
that , 


1// 1 + R2 


so 


F (R)  = 


dR’ 


° / T2 

/ 1 + R' 


= sinh  ^ R 


(BA) 


hence  a = a sinh  R ; R = sinh(a/a  ) 
o o 

Frim  Equations  (B4)  and  (B5)  we  thus  have 


F’(R)  - 


/ 1 + r" 


cosh(a/a  ) 

O 


Finally,  substituting  Equation  (B6)  in  Equation  (B3)  the  combined  plane- 
spherical  wa^e  form  of  Burgers’  equation  expressed  in  terms  of  'stretched' 


coordinates  becomes, 


— _ (l/N)  W — 

3o  U/  ' 8t 


T cosh(o/o  ) 
o o 


1 

1 


For  a monofrequency  source  (i.e.  N = 1)  the  solution  of  Equation  (B7)  can 
be  expressed  via  the  Fubini,  Fay,"^  ^ and  asymptotic  far-field"^  ’ 

approximations  as. 


00  2J  (na) 


W(a,T) 


sin(nT)  , 


o < a < 1 (B8a) 


„ 00  cosh(o/a  ) sin(nl) 

_L  i 

P ^ “ 1 

o n=l  sinh{nr  (1  + a)  cosh(o/a  )} 
o o 


, tt/2  < o < o (B8b) 
— — s 


9 MA> 

j y-— y (exp[-(oo/ro)  sinh(o/ao)]}sin(T)  , 0g  < 0 < °°(B8c) 
o 


where  A = (2/r^)  coshCa^/a^) 

and  the  'scaled'  shock-termination  distance  a is  obtained  by  solving 


transcendental  equation,  previously  expressed  by  the  author  in  slightly 
different  form  for  plane  or  spherical  waves  as, 


r i 0 i 

0 0 

cosh(o  /a  ) = 

c o 

f ' 

x r 

n n 

o o 

o U 

u u 

L i 

By  varying  Oq  we  can  compute  °s/°0  f°r  particular  values  of  aQr0  » 
evaluate  A via  Equation  (B9)  and  thus  obtain  the  total  finite-amplitude 
absorption  loss  EXDB^  for  Equation  (B8c)  in  the  form. 


EXDB 

00 


20  L°SlO 


2 

A Io(A)  ' * 


Values  of  EXDB  obtained  in  the  manner  are  found  to  be  in  good  agreement 

00 

with  similar  values  calculated  with  the  aid  of  Equation  (31b). 
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